Electronic spin drift in graphene field effect transistors 
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We studied the drift of electron spins under an applied DC electric field in single layer graphene 
spin valves in a field effect transport geometry at room temperature. In the metallic conduction 
regime (n ~ 3.5 x 10 16 m -2 ), for DC fields of about ±70 kV/m applied between the spin injector 
and spin detector, the spin valve signals are increased/decreased, depending on the direction of the 
DC field and the carrier type, by as much as ±50%. Sign reversal of the drift effect is observed when 
switching from hole to electron conduction. In the vicinity of the Dirac neutrality point the drift 
effect is strongly suppressed. The experiments are in quantitative agreement with a drift-diffusion 
model of spin transport. 
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Graphene, a two-dimensional hexagonal carbon lat- 
tice, was isolated from graphite and became available 
for electric transport type measurements in 2004 [l], Q • 
Graphene is unique in its nature by allowing the exper- 
imentalist to shift the Fermi level trough the zero-gap 
energy bands and continuously tune the carrier density 
from holes to electrons via the Dirac neutrality point. 
Recent studies also addressed the spin degree of free- 
dom, resulting in successful injection, transport and de- 
tection of spins in single and multiple layers of graphene 
@) BL H, B S H] • The spin relaxation time r s f measured 
about 150 ps (as extracted from Hanle type spin pre- 
cession measurements) , which combined with a diffusion 
constant of D ~ 0.02 m 2 /s yielded a considerable room 
temperature spin diffusion length X s f = ^JDt s j of 1.5 to 
2 /mi. 

Until now, spin transport in graphene was studied in 
the low bias regime where the transport is described by 
diffusion only. In this paper we study the spin trans- 
port in the presence of DC electric fields. Due to the 
large carrier mobilities in graphene [9j], the DC fields give 
rise to spin drift effects which are comparable in mag- 
nitude to the spin diffusion in graphene and, therefore, 
become accessible experimentally. In order to study the 
spin drift effect, we have built lateral spin valve devices 
as illustrated in Fig. Q] The graphene layer is contacted 
by a series of Co electrodes (F1-F5) used for the injec- 
tion/detection of spins and the application of the DC 
electric fields. The dimensions of the contacts and their 
spacings given in the figure caption are the ones of the 
sample discussed in this paper and are characteristic for 
the spin valves we have fabricated. The graphene flakes 
we use are typically between 0.2 and 4 fim wide and up to 
50 /jm long. Besides the AFM thickness characterization 
of every flake, on several samples we have conducted Ra- 
man spectroscopy measurements as well [ll| , in order to 
select the single layer flakes. For a detailed description of 
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the device fabrication see Ref. [||. The contact resistances 
are between 2.5 kfi and 12 kil. A standard low-frequency 
AC lock-in technique with a current of 4 /iA RMS was 
used for measuring the spin signals. The measurements 
were performed at room temperature and in vacuum in 
order to minimize the shifts of the Dirac point due to 
charging effects of adsorbates on the graphene layer. 

An AC spin imbalance is created in graphene by in- 
jection/extraction of spin currents by the ferromagnetic 
electrodes F4 and F5 through an aluminium-oxide tunnel 
barrier. Due to the opaque C0/AI2O3 /graphene inter- 
face, the spin transport is not disturbed by the central 
contact F3. The spin imbalance diffuses symmetrically, 
with respect to the injection points F4 and F5, through 
the device giving rise to an AC voltage drop between the 
ferromagnetic detector electrodes F2 and Fl. In this non- 
local measurement geometry, the spin transport is thus 
completely separated from the charge transport and is 
solely described by the diffusion constant D and the spin 
relaxation time r s f . The non-local voltage drop between 
F2 and Fl is measured as a function of a magnetic field 
applied along the easy axis of the Co electrodes (y direc- 
tion). This results in typical spin valve signals as shown 
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FIG. 1: Schematic drawing of the spin drift experiment. The 
Co contact widths are 350 nm (for Fl), 90 nm (F2), 250 nm 
(F3), 90 nm (F4), 350 nm (F5) and the distances between 
them are (from left to right) 2 ^m, 0.6 fim, 1.3 /im and 2 /jm. 
The 0.9 x 7 /im 2 sized graphene flake is covered by an alu- 
minum oxide layer. A (gate) voltage, applied between the Au 
contacted Si(n ++ ) substrate and the grounded electrode F5, 
allows for the control of the charge density in graphene. The 
layer thicknesses are given in the brackets. 
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in Fig. [2] a). The variable width of the Co electrodes 
yields a difference in their coercive fields (20 to 50 mT in 
the present case), allowing for separate switching in the 
external magnetic field. Contacts with similar widths, 
however, can sometimes interchange in their switching 
order depending on, e.g., the domain wall nucleation, on 
the sweeping direction of the field (positive to negative 
values and vice versa). Such an effect is observed in the 
measurement shown in Fig. [2] a) where the switching or- 
der of F4 and F5 is interchanged on the negative versus 
positive field's side. The identification of the switches 
could nevertheless be done by examining the spin valve 
signals obtained in different regions of the same sam- 
ple by choosing different pairs of injector/detector elec- 
trodes. 

Each discrete resistance level seen in Fig. [5] a) 
corresponds to a combination of magnetically paral- 
lel/antiparallel electrodes. The step heights carry infor- 
mation about the importance of that specific electrode's 
contribution to the spin injection/detection process. The 
highest change in signal is observed when the inner elec- 
trodes (F2,F4) change their relative parallel/ antiparallel 
magnetization state. The further away an electrode is, 
the smaller its corresponding resistance step is. The spins 
injected at F5 must diffuse longer (and have thus higher 
chance to relax) before they are detected by the detectors 
(F2,F1) than the spins injected at F4. The presence of 
four switches in our measurements proves the spin trans- 
port through the graphene channel under all magnetic 
electrodes, over the full ~ 6 fim length of the device. 

Spin drift effects in all-electrical n-GaAs based devices 
were measured by Lou et al. [10t | . where the conduc- 
tion was limited to electrons. Unlike in semiconduc- 
tors, in graphene, due to its band structure, electro- 
static gating allows for switching from hole to electron 
conduction while keeping the carrier mobility, diffusion 
constant, Fermi velocity, electric conductivity and other 
parameters approximately unchanged. The gating ef- 
fect is reflected in the graphene resistivity as plotted in 
Fig. [2Jd) • The position of the Dirac neutrality point sepa- 
rating hole and electron conduction regimes, as identified 
by the minimum in conductance (l2l . Il3j , shows a small 
hysteresis (+ 9 V to + fl V) for the two sweep direc- 
tions. The measurement also reveals a carrier mobility 
of fi =0.25 m 2 /Vs. 

Applying a DC electric field to generate carrier and 
spin drift in addition to the diffusion would in princi- 
ple be possible by passing a DC current from F2 to F4. 
However, a strong DC bias can dramatically influence the 
spin injection efficiency of the electrodes [14j, [l5| mak- 
ing it impossible to separate the drift effect contribution. 
We have conducted experiments on graphene spin valve 
devices where both AC spin injector electrodes are DC 
biased and we measured changes of the spin valve effect 
up to an order of magnitude compared to the non-biased 
measurements [1 61 ] . In order to avoid the manipulation of 
the spin injection efficiency, we have introduced the elec- 
trode F3 (see FigfTJ). The DC field is applied between the 
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FIG. 2: a) Non-local spin valve measurement of the device 
presented in Fig. [1] The small arrows represent the magnetic 
orientation of electrodes Fl, F2, F4 and F5. The asymmetric 
spin valve behavior is due to the switching order of F4 and F5 
which depends on external field sweep direction (indicated by 
the large horizontal arrows). Rsv represents the spin signal; 
see text for details, b) 4-point measurement of the graphene 
resistivity between contacts F3 and F4 versus the gate voltage. 



electrodes F3 and F5 by sending a constant DC current 
between them, with F5 being grounded. We could not 
avoid sharing the drain of the AC and DC circuits, the 
injector F5 on the far right of the device is effectively DC 
biased. However, the effect of changing its injection effi- 
ciency with the DC current bias is minimized as discussed 
later. 

Experimentally, the drift will manifest itself as an 
asymmetric distribution of the spin imbalance on the two 
sides of the spin injection point F4 (and F5), in addition 
to the symmetric diffusion process. Between F3 and F5 
thus, where the field E is active, a favored direction of 
the spin transport will be defined by the electric field's 
polarity and the carrier type (electrons or holes). Such 
a device with fixed-distance injector/detector electrodes 
will therefore show a variable spin signal dependent on 
the DC electric field. 

In Fig. [3] we present room temperature spin valve mea- 
surements for several DC currents (applied between F3 
and F5) at three gate voltages: -20 V (hole conduction), 
+9 V (the Dirac neutrality point) and +38 V (electron 
conduction). The upper limit of the electric fields was 
set by the current densities allowed to pas s through the 
tunnel barriers without damaging them[l8j. Before each 
measurement, the magnetic field is ramped up to 300 mT 
to saturate the magnetization of the contact electrodes. 
The DC currents used were ±40, ±20 and /iA and their 
corresponding electric fields were ±68 kV/m, ±34 kV/m 
for V s =-20 and +38 V and ±206 kV/m, ±103 kV/m for 
V g =+9 V. Similar measurements, done on a set of sam- 
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pies with different geometry and/or contact resistances, 
support the data we present and discuss here. 

The change in the injection efficiency of the electrode 
F5 (common ground for the DC and AC currents) for dif- 
ferent DC fields appears in the spin valve measurements 
as an increased/decreased contribution of this electrode 
to the device non-local resistance. This can be seen in 
Fig. [3] where, for the measurements performed at +38 V 
gate voltage, the resistance step at approximately -30 mT 
(corresponding to the magnetic switch of electrode F5) 
does show a dependence on the DC electric field. How- 
ever, this electrode is positioned far from the AC detec- 
tion circuit, and we can minimize its effect by defining 
the "spin signal" Rsv as the resistance step between the 
parallel and antiparallcl magnetic state of the two cen- 
tral electrodes F2 and F4 (see Fig. [5]). Therefore, the DC 
field dependence of Rsv indicates the presence of carrier 
drift as explained below. 

At -20 V gate voltage (left column of graphs in Fig. [3]) , 
we are in the high density hole conduction regime with 
a carrier concentration of rih — 3.5 x 10 16 m -2 yielding 
a graphene resistivity R sq uare — 1-5 kf2. The non-local 
spin signal Rsv reads 1 SI when no DC field is present 
(graph c). Applying a DC field of +34 kV/m (graph d) 
and +68 kV/m (graph e) decreases the spin valve signal 
to 0.7 ft and 0.45 ft respectively, indicating spin drift 
in the region F3-F5 that counteracts the symmetric spin 
diffusion. This transport is labeled as "upstream" since, 
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FIG. 3: The drift effect in a graphene spin valve. The gate 
voltage V 9 varies from left to right and it is indicated for each 
column. The DC electric current varies from top to down as: 
(a) I DC = -40 fiA, (b) -20 uA, (c) fiA, (d)+20 fj,A, (c) 
+40 fiA. For clarity only the negative side of the spin valve 
measurements are plotted. 



in order to detect a spin valve signal at electrodes F2 
and Fl, the spins need to travel against the action of the 
carrier drift. Now, if the DC electric fields' polarity is 
reversed, i.e. we apply a field of -34 kV/m (graph b) and 
-68 kV/m (graph a), the resulting drift facilitates spin 
transport towards the detectors ("downstream" trans- 
port) and the measured spin valve signal increases to 1.3 
and 1.5 17 respectively. 

In the case of electron conduction, the drift effect re- 
verses, because the electric field - charge carrier interac- 
tion is opposite, as seen in the set of curves measured 
at V 5 =+38 V (right column in Fig J3]). This voltage 
was selected to give approximately the same parameters 
(n e , Riocai) on the electron conduction side of the Dirac 
neutrality point as for the hole conduction (V g =-20 V). 
The spin signal scales with the field amplitude and direc- 
tion as expected, reading 0.6, 0.9, 1.2, 1.5 and 1.7 ft for 
-68, -34, 0, +34 and +68 kV/m respectively. 

At a gate voltage that locates the Fermi level close to 
the Dirac neutrality point (R sqU are — 4.5 kfi) the conduc- 
tion happens neither via electrons, nor via holes, which 
means that no spin drift should take place. Examining 
the measurements at V g =+9 V (central column of graphs 
in Fig. [3]), we observe a slight decreasing tendency of the 
spin valve signal for increasingly positive DC fields, simi- 
lar to the drift in case of hole conduction. This behavior 
can be explained by the interplay of the following two 
effects. 

First, under the high DC current a certain shift of the 
effective gate voltage is expected. Considering the 20 kfi 
resistance of the structure (graphene + contact F5), this 
shift equals to ±0.8 V for currents of ±40 /iA. Unlike in 
the case of high carrier densities, around the neutrality 
point such a shift can critically change the carrier density 
in graphene and lead to spin drift. Second, due to a shift 
of the Dirac point in time under a constant gate voltage 
stress, it is difficult to set and keep a gate voltage just 
at the right position to assure neutrality. In our case, it 
seems that we were slightly in the hole conduction regime 
when measuring at V g =+9 V. 

For a quantitative interpretation of the data we have 
extracted the spin signal values from all the measure- 
ments of Fig. [3] and plotted them in Fig. 0] against the 
DC field for the three gate voltages. For the theoreti- 
cal description, we adopt the drift-diffusion model intro- 
duced by Yu and Flatte [3, QjJ. In steady-state, the 
drift-diffusion equation for the spin imbalance n s reads: 



DV 2 n s + fiEWn a 



0. 



(1) 



where D is the diffusion constant, r s / is the spin relax- 
ation time, /i is the carrier mobility and E is the DC 
electric field. The term [iE represents the drift veloc- 
ity Vd and its magnitude compared to the Fermi velocity 
Dp = 10 6 m/s is important for the contribution of the 
spin drift to the spin transport. Similar to the spin dif- 
fusion length X s f = yj ' Dr s f, we define a spin drift length 
A f ; = D/vd- The symmetric diffusion and asymmetric 



4 



drift effects add up to form a spin transport character- 
ized by a pair of len gth scales, named upstream and down- 
stream lengths A± fl71 ]: 



f 



11 /ll 1 



(2) 



The general solution to the spin imbalance equation [T] in 
the direction x parallel to the spin transport is 

n s (x) = Aexp(+-^) + £exp(-^-), (3) 

A+ A — 

where A and B are determined by the boundary condi- 
tions. Since we take the electrode F4 as the spin source, 
only the decaying solution with respect to the injection 
point need to be considered. Therefore, we can write the 
spin signal (Rsv) dependence versus the applied DC field 
in the form: 



Rsv = #oexp(- — ), 



(4) 



where Ro is the spin signal that would be measured at the 
injection point, A± carries the electric field dependence 
(Eq. \2§ and L ~ 1.5/im is the distance between F3 and 
F4. 

In the following we compare the drift-diffusion equa- 
tion for the spin signal (Eq. 0]) with the spin drift mea- 
surements in the high carrier density regime. From the 
conductivity measurements presented in Fig. 2b we ex- 
tract /i ~ 0.25 m 2 /Vs using a — nefi, where the carrier 
concentration n is calculated from the sample's response 
to gate voltage (fig. 2b), using the method described in 
Refs. I2GJ. This gives a drift velocity of v d ~ ±1.7 • 10 4 
m/s for E ~ ±68 kV/m. The diffusion constant of 
D ~ 0.02 m 2 /s is obtained using the Einstein relation 
(7 = Ne 2 D, where the density of states N is determined 
according to Ref. H Using Eq. [2] and a spin relaxation 
length X s f ~ 2^m [H, we calculate the magnitude of 
the drift effect manifesting in the up/downstream lengths 
as A + = 1.3/im, A_ = 3.0/im for E =±34 kV/m and 
A+ = 0.9^m, A_ = 4.4^m for E =±68 kV/m. Note, 
that at E =±68 kV/m the characteristic length scale 
of the drift-diffusion spin transport is about two times 
increased/decreased compared to the diffusion process. 
The value of Rq we estimate from the spin valve mea- 
surements in zero electric field when X± — A s / and Rq is 
given by i?o = Rsv exp (L/A s /). By taking the length of 
the drift region L = 1.5/xm, A s / = 2/um and the measured 



Rsv at E =0 we can determine Rq. Now, we are able to 
calculate the drift effect according to Eq. 4, all parame- 
ters being known. The theoretical curves (solid lines in 
Fig. [J) are in excellent agreement with the experimen- 
tal data. There are no free parameters in the calculation 
and, in fact, these measurements could be used to extract 
A s /, if it would be unknown. 

At the Dirac point, due to the symmetry, no drift effect 
is expected. However, in our measurements we are prob- 
ably slightly away from the Dirac point. Here we cannot 
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FIG. 4: Spin valve signal for high carrier densities (circles), 
and near the Dirac neutrality point (stars). The lines repre- 
sent a calculation (with no free parameters) of the drift effect 
based on the drift-diffusion equation. 



apply our above analysis. It remains an open theoreti- 
cal question to what the drift velocity should be at the 
charge neutrality point. 

In conclusion, we studied spin transport manipula- 
tion in graphene, by carrier drift under the action of 
a DC electric field. For high charge carrier densities 
(n ~ 3.5 x 10 16 m~ 2 ) depending on the direction of the 
applied DC field and the nature of the carriers we were 
able to modify the effective spin relaxation length by fac- 
tor of 4.8 which resulted in a modulation of the spin valve 
signal of about 300%. The spin-drift measurements are 
described well by a drift-diffusion model. The control 
over the drift velocity we demonstrated opens new pos- 
sibilities for exploring spin related phenomena in other 
types of graphene devices, such as graphene p-n junc- 
tions. 
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